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Abstract 

The corrections for two-pion correlations due to electromagnetic final-state inter- 
actions at high secondary multiplicities are investigated. The analysis is performed 
by solving the Schrodinger equation with a potential which is dictated by the multi- 
particle environment. Two different post-freeze-out scenarios are examined. First, for 
a uniformly spread environment of secondary particles, a screened Coulomb potential 
is exploited. It is shown that the presence of a static and uniform post-freeze-out 
medium results in a noticeable deviation from the standard Gamov factor. However, 
after going to a more realistic model of an expanding pion system, this conclusion 
changes drastically. We argue that the density of the secondary pions n n (t, R), where 
R is a distance from the fireball, is bounded from above by n n (t,R) < const/i? 2 for 
all times t. Then, a two-particle scalar potential which is found as a solution of the 
Maxwell equation for non-uniform medium replaces the screened one. Even this upper 
limit does not result in an essential deviation from the Gamov correction. 
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1 Introduction 



Two-particle correlations provide information about the space-time structure and dynamics 
of the emitting source [[L]]. Considering the correlations that occur in relativistic heavy- 
ion collisions one usually assumes that: (i) the particles are emitted independently (or the 
source is completely chaotic), and (ii) finite multiplicity corrections can be neglected. Then 
the correlations reflect a) the effects from symmetrization (antisymmetrization) of the wave 
function and b) the effects that are generated by the final-state interactions of the detected 
particles between themselves and with the source. At first sight one can regard the final- 
state interactions (FSI) as a contamination of the 'pure' particle correlations. It should be 
pointed out, however, that the FSI depend on the structure of the emitting source and thus 
provide information about the source dynamics as well. In fact, this was proved by intensive 
investigations during the last twenty years of two-particle and source-particle FSI [[!], 0] (for 
recent publications see, for example, on two-particle FSI, ||, |J and references therein; on 
source-particle FSI, [|5], |5| and references therein). 

Actually, the former approaches accounting for FSI deal with the secondaries in empty 
post-freeze-out space. Meanwhile, in recent SPS experiments, for instance Pb + Pb at 
160 GeV x A, some 800-900 secondary charged pions are created which form obviously 
a plasma-like post freeze-out medium. Therefore, one might expect that the FSI of two 
separate pions, at this collision energy, and of course at the energies of the forthcoming 
RHIC and LHC colliders, would be strongly influenced by the environment formed by other 
particles. The goal of the present paper is to estimate how large are the consequences on the 
Coulomb final-state interactions due to the presence of a large number of secondary charged 
particles. To be as independent as possible from source models, we chose the estimation of 
the Gamov factor as a standard quantity, which serves as a measure of the Coulomb FSI. 

The fundamental observable for intensity interferometry in hadron physics is the relative 
momentum spectrum of identical particles. For two like-sign pions, the modifications to the 
spectrum by the final-state Coulomb interaction result in a correction that has typically been 
considered as tractable with high accuracy. This assumption was based on the significantly 
different length scales between strong (oc l/m n ) and Coulomb (oc l/m n a) interactions J7|. |J 
(here m n is the pion mass and a represents the fine structure constant). The correction may 
then be treated on the basis of the Schrodinger equation, resulting in the well-known Gamov 
factor G(q) 

G(|q|)=|^(r = 0), e27rj? _ i; 

where T/> q (r) is the two-particle wave function, rj = am n /\q\, and q is the relative momentum 
of the particles. 

The nominal quantity expressing the correlation function in terms of experimental dis- 
tributions U is 

where Pi (k) = Ed 3 N/d 3 k and P 2 (ki,k 2 ) = E x E 2 d 6 N/(d 3 k 1 d 3 k 2 ) are single- and two- 
particle cross-sections. For point-like emitters it can be expressed (due to the factorization 
of the corresponding matrix element) in terms of a product of the Gamov factor to the model 
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correlations [0, [J: 

C(k 1 ,k 2 ) = G(|k 1 -k 2 |)C modcl (k 1 ,k 2 ). (3) 

It should be mentioned that, in the centre of mass of the pair the relative pion momentum 
|q| = 2|k|, with k = ki = — k 2 , coincides with invariant relative momentum qw = \{k\ + 
A; 2 ) 2 — 4m 2 ] 1 / 2 , where k\ and /c 2 are the pion four-momenta in an arbitrary frame. 

We consider corrections due to the Coulomb final-state interactions in a common scheme 
in which we do not take into account the finite-size of the fireball (for a consistent treatment 
of the source finite-size effects, see for example ||). This means that we calculate the 
correction factor in accordance with the formula G C orr(|q|) =| V^ 1 " = 0) | 2 , where now ^ q (r) 
is obtained from a numerical solution of the Schrodinger equation with the two-particle 
Coulomb potential distorted by a multipion environment. Since, as shown in ||, the finite 
size of the emission source softens the manifestation of the FSI, the 'Gamov factor' tends 
to overestimate the FSI effects and therefore it is the most sensitive quantity for deviations 
from the standard two-particle Coulomb interaction. 



2 Static multiparticle environment 

In the high multiplicity case, when a post-freeze-out multipion environment cannot be ne- 
glected, the relation between the two-particle electromagnetic potential 0(r) and the local 
charge density is given by 

V 2 0(r) = -Ane(n {+) - n H ) , (4) 

where e = yfa is the elementary charge; the density of charged pions rv^' is related to that 
of neutral pions via a Boltzmann factor: 

n (±) = n (o) exp ; (5) 

The density of 7r°-mesons at the freeze-out temperature Tf coincides with the equilibrium 
density of charged pions in the absence of Coulomb interactions (we consider symmetrical 
nuclear matter). In the limit e0 Tf, Eq.([|) can be rewritten as 

n (±) = - (0) (l T |) (6) 

(this requires that the pions are not closer than ~ 10~ 2 fm to one another at Tf « 200 MeV), 
so that 

V 2 0(r) = — (2n<°>Mr) . (7) 
The solution of this equation is well known and given by a screened Coulomb potential 

e" 



-r/i? scr 



07r±( r ) = ±e- 

where 



-^scr V 3 
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a ' \ 7FT (9) 
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with the total pion density = 3n^°\ Thus, for like-sign pions the potential energy reads 

U 7T7T (r) = -exp(--^-) . (10) 

To evaluate the correction factor we use the screened Coulomb potential fllPf ) to solve 
the Schrodinger equation numerically, for the two choises: 1) = 0.25 fm -3 and 2) = 
0.03 fm~ 3 . (From now on we shall quote these two cases as "Li/C" and "5"P5" freeze-out 
conditions, respectively.) Taking Tf = 190 MeV, for example, we obtain from Eq. (g) 
screening radii i? scr ~ 7.9 fm (LHC), R scr ~ 22.4 fm (SPS), respectively. For completeness 
of illustration we also consider the intermediate case _R scr ~ 19.3 fm. 

The results of these calculations are plotted in Fig. 1 together with the standard Gamov 
factor. We see that a substantial correction to the standard Gamov factor would be required 
even for the existing experimental data if one adopts the idealized picture of a uniform 
post-freeze-out density of the environment. 

It is interesting to point out that the correction factor G COI (q), which was evaluated for 
the same screened potential using the quasi-classical approximation |J, does not make a big 
difference with the correction factor in Fig. 1. 

Actually, what we really learned from the consideration of this idealized scenario is that 
a cancellation of the Coulomb potential tail (by screening) results in an increase of the 
correction factor in the region of small relative momentum (q < 50 MeV), as can be seen 
from Fig. 1. Thus the long-distance behaviour of the potential is responsible for the dramatic 
deviation of the correction factor from the standard Gamov factor at small relative momenta. 
On the other hand, the tail of the two-particle potential gives the main contribution to 
interactions when the particles are at large distances from one another, where in turn the 
density of secondary particles is small in the realistic picture. To keep this qualitative 
speculation as a thread we next turn to a more realistic calculation, explicitly incorporating 
expansion. 



3 Expansion scenario 

In the previous scenario the whole position space was filled by particles with a constant 
density; this is certainly an unrealistic approximation (idealization). In order to take into 
account post-freeze-out expansion of the pion system, we parametrise the pion density as 

n{R) = nt£, (11) 

where rif is the freeze-out pion density and Rf is the freeze-out radius. Indeed, the spatial 
volume of the expanding pion system in the solid angle Q increases as AV = Q ■ R 2 ■ AR, 
where R is the distance from the centre of the fireball and AR is the thickness of the layer, 
which we keep constant. Then, if the number of particles AN in this volume is constant, 
the density reads: n(R) = AN/(Q ■ R 2 ■ AR) = const/i? 2 . Thus, the model (|TT|) implies that 
all particles have the same modulus of radial velocity. As we shall see further, such a model 
still means an overestimation of the particle density. 
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To support our assumption (|TT|) let us consider a classical pion phase-space distribution. 
After freeze-out it satisfies the collisionless Boltzmann kinetic equation 

^^ + v.V/(x )P ) = 0, (12) 

where v = p/po is the velocity of the particle and po = ^(p) = ym^ + p 2 . We look for the 
expansion solution of this equation, which can be fixed by asymptotic condition, for instance 
lim^oo f(t, x = 0; p) = 0. A solution of this type can be written in the form 

/(f,R,p) = /o(R-v*,p) , (13) 

where / (R, p) is the initial distribution (t = 0). For the sake of simplicity, we take an 
isotropic initial distribution in position and momentum spaces: just before freeze-out the 
particles were distributed (a) in accordance with Boltzmann's law in momentum space, 
and (b) in accordance with a Gaussian distribution in position space. The classical kinetic 
equation is quite sufficient to describe the collective (!) behaviour of the pion system after 
freeze-out. Hence, we assume that the system, at time t = 0, occupies the phase space 
according to the distribution function 

/ (R,p) = n (R)<7o(p) , (14) 

where 

with / d 3 Rn (R) = N n , and 



M R ) = 7Z ^lv 2 ex P (-^2) ( 15 ) 



WP> = /m x exp - V ] (10) 



<T { K 2 



with J[d 3 p/ (2ti) 3 } g (p) = 1. Here N n is the total number of pions, Tf and R{ are the 
temperature and the mean radius of the system at time t = (freeze-out), respectively, m n 
is the pion mass and K2 is a Bessel function of imaginary argument. 

The spatial distribution of the particles at time t is determined by integrating the distri- 
bution function (|H|) over the momentum variable 

n(t,R)= f -fP-n (R--?-t)g (p) . (17) 



(2tt)3 u V w(p) 

Because of the spherical symmetry, it is reasonable to look at the radial density of pions 

n sph (t, R) = AnR 2 n(t, R) . (18) 

This quantity may be treated as the number of pions in the shell with unit thickness at time 
t and at a distance R = |R| from the fireball centre. Hence, n sp h(t,R) is a one-dimensional 
spatial distribution function and, evidently, the area under this curve at any time is equal 
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to the particle number N n , because of its normalization J °° dRn sph (t, R) = N w . We evaluate 
this function in accordance with Eq. ( |17|) for SPS freeze-out conditions: rif = 0.03 fm -1 ^ 3 , 
Tf = 190 MeV and i?f = 7.1 fm. The results of this calculation at different times t are 
given in Fig. 2. It shows that the spherical distribution is always almost Gaussian-like 
and the velocity of the distribution maximum is very close to the velocity of light. The 
horizontal line (see Fig. 2) denotes a constant spherical density 47rR 2 n(R) = const. This 
line is nothing more than the 3-dimensional spatial density n(R) = const'/i? 2 (for the SPS 
freeze-out conditions const' ~ 85/47r). It follows that, at any time t, for the post-freeze-out 
particle density n(t, R) the upper limit is valid 

n(t,R)<n^, (19) 

where rif = max{n(t = 0,R)} and the equality is reached for an expanding system where 
the radial particle velocities are equal. Because of momentum dispersion in the post-freeze- 
out pion system, the pion density (|TT|) is an essential overestimation of the real density, 
which is formed by the comoving multipion environment. Hence, adopting the stationary 
density dependence (|ll|) we overestimate the influence of the environment, and consequently 
overestimate a distortion of the Coulomb FSI in the expansion scenario. On the other 
hand, this approach provides us with a possibility to consider a stationary post-freeze-out 
environment even in the frame of the non-stationary expansion scenario (see Fig. 2). 

Adopting the parametrization (JTT|), we get Eq. (^), where the pion density n n now depends 
on R: 

V 2 0(r) = f^n(R) 0(r) , (20) 

where we put « n n /3 as before (rif = n n ). Inserting the pair centre of mass distance 
R together with the distance r between two detected particles, we have, in the classical 
approach: 

Rm Rf + v cm ■ t , r&Vrd-t, (21) 

where v cm is the velocity of the two-particle centre of mass in the fireball rest frame and 
f re i is the relative velocity of the particles (t> re i = q/m, t = is fixed on the freeze-out 
hyper-surface). Eliminating the time from the approximate equalities ( pl|) one has 



R = R { + . r , (22) 

Vrel 

which means that we parametrize the time evolution by the mean distance r between two 
detected particles. 

We can thus rewrite the dependence of the pion density on r as 

»W) = (-) ! n^. P3) 

\v cm J (r + r) 2 
where we define the dynamical freeze-out radius 

r = Rf ^1 . (24) 



v, 



cm 
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It is time to recapitulate what has been done. First, we eliminate the time dependence of 
the particle density, bounding it from above by a stationary non-uniform spatial distribution, 
in accordance with inequality (Jl9|). Then, we connect the mean distance R of the pair c.m.s. 
from the fireball and the mean distance r between pions, thus eliminating again the time 
dependence. This means that all time and R dependences are now parametrized by the 
variable r. If the relative velocity of the separate pions is small, then the partial time 
derivative = v re \ can be neglected. Hence, the problem is approximately reduced to a 
stationary one. Solving Eq. (|20|) with particle density from fl23|), one obtains the two-particle 
potential energy U(r) = e0(r), which will then be exploited in the stationary Schrodinger 
equation to find a wave function. This will be the further strategy of our estimations. So, 
the problem of a time-dependent (expansion scenario) is reduced to a stationary one by the 
price of somewhat overestimating the Coulomb corrections. 

Equation ( p0|) may be rewritten as 

V72^ _ c 2 (g) ^ (25) 



(r + r) 2 ' 

where we note explicitly that when the particles are separated by large distances r the density 
of the multiparticle environment goes down in accordance with (|23[) . The quantity c 2 (q) is 
defined in the following way: 

^ = MaRjr^v^ 



cm 



Fixing the screening radius at freeze-out 

<27) 

we have 

c( ? ) = • (28) 

Rscr "cm 

In spherical coordinates, Eq. (|25| ) can be rewritten as 



#m + im.- *<sL« T) = . (29) 



dr 2 r dr (r + r) 2 
This equation may be solved by 



(30) 



where 0(r) satisfies the boundary condition: 0(r) — > e/r when r — ■> 0. The exponent b is 
a solution of the quadratic equation resulting from the substitution of the ansatz (^) into 
Eq. ( p9| ) and takes the form (assuming proper asymptotic behaviour of the potential 0): 



%) = -2 + 2V 1 + 4c2 (<?)- ( 31 ) 

One has 6^0 when ri{ — > and hence the potential 0(r) transforms into the Coulomb 
potential in this case. It is shown in Figs. 3a and 3b for LHC and SPS freeze-out conditions, 
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respectively. In the interval of interest, b increases with increasing relative velocity, and so 
do deviations from a pure Coulomb field. This intriguing behaviour is directly related to the 
"Hubble-like" expansion implied by Eq. fl2"2"|): it becomes clear if we remember that the pion 
density decreases (hence R scr — > oo) with increasing distance R. The expansion thus results 
in modifications to the Coulomb potential that are of power-law, not exponential, form, in 
contrast to the static result given by Eq. (H). 

The connection of the potential (^) and screened one can be found in the following way. 
One may treat the corrected potential obtained in Eq. (|30|) as an effective charge distribution 



e e g = e — - , (32) 
\r + rj 

which we are going to average. Equation (^) can be rewritten as (V 2 — K 2 )<fr(r) = 0, with 
k = c(q)/(r + r). This is equivalent to an r-dependent screening radius, i.e. 

r + r 

Rscr(r) = —r^- • (33 
c{q) 



As shown in Figs. 3, the deviation of the potential (|30|) from the pure Coulomb form in the 
region of small relative pion momentum q < 30 MeV, is small (b — > 0). Since, k —>■ when 
v re [ <C 1 (see Eq. (p8|)), we first ignore the r dependence of k to obtain the solution for the 
electromagnetic irn potential in the form of Eq. (§), then substitute the r-dependent k into 
Eq. (|), to find that 

ae -c(q)r/{r+r) 

= (34) 

r 

no longer decreases exponentially with distance when r is large enough r 3> r. Instead, the 
numerator on the r.h.s. of Eq. (0) represents the averaged charge distribution (^) squared 
(we are now considering the potential energy U n7r rather than the electric potential 0, hence 
the extra factor of charge leading the a). 

It is clear from Eq. ( pgp that if v re \/v cm 1 (i?f and R{ cr are of the same order for high 
multiplicities) the renormalized constant a e s = a exp [— c(q)] is close to the bare value of a. 
Moreover, the same qualitative result comes from the r-behaviour of the screening radius 
( |3~3"D when it approaches the asymptotic value i? scr = oo (Coulomb law) with increasing r. 
The quantity c(q) increases with relative pion momentum, leading to larger deviations from 
the Coulomb potential and agreement with the features of the potential (|3~0"D (see discussion 
after Eq. ©). 

3.1 Evaluations 

The numerical evaluations of the correction factor (Gamov factor) should be provided by 
the solution of the Schrodinger equation with the potential 

t/cxpan(r) = - (^—) b , (35) 



where the exponent b is momentum-dependent, in accordance with fl3"T|). Then, one can 
construct the correction factor G C or(|q|) =| ^q( r — 0) | 2 . To compare it with the standard 
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Gamov factor G we calculate the ratio Gq/G cot for the LHC freeze-out conditions. The 
results are depicted in Fig. 4 as dotted curves for three different values of the pair mean 
momentum P cm . 

1/3 

We now correct the two-particle potential (|35| ) for small distances r < a = n { , where 
rif is the freeze-out pion density, i.e. in the region where the distance r between two pions 
is smaller than the mean distance between the particles in the gas. Indeed, there is no 
'screening' effect for this distance, and hence the potential should be the Coulomb one. 
There is no such a problem for the statically screened potential U{r) = aexp (— r/R SCT )/r, 
because the screening radius obeys the condition a -C R SCT (by definition in the sphere 
47r_Rg Cr /3 the number of particles is much larger than 1); for the region r < a <C _R scr , the 
potential U(r) automatically transforms into the Coulomb one, U(r) = a/r. In the expansion 
scenario we have in addition to R SCT another scale parameter, f = Rf (v ie i/v cm ), where R{ is 
the freeze-out size of the system (i? f = 7 — 10 fm for SPS and LHC freeze-out conditions). 
By construction the potential U expaun (r) approaches the Coulomb one if (i) r C r (see (pop). 
On the other hand, the scale parameter r may be smaller than the mean distance between 
particles in the gas (ii) r < a. In Fig. 5 we depict r for different mean momenta of the 
pair: P cm = 50, 150, 450 MeV/c and for the mean distance between particles at freeze-out 
a = n { (horizontal lines). For instance, as is seen in Fig. 5, for P cm = 450 MeV/c and 
relative momentum q < 60 MeV/c, we have r < a for SPS freeze-out conditions. Hence, 
combining (i) and (ii) we get that the asymptotic regime, i.e. the Coulomb potential, can be 
achieved only for r <C a. But we know that for separation distances which obey r < a, all 
distortions of the Coulomb potential vanish. Thus, the behaviour of the potential ( |35"D must 
be corrected, since it should be a Coulomb one at distances that are smaller than the mean 
distance between particles. To improve the behaviour of the potential for distances smaller 
than a, we use a smooth 'potential switcher' 



which is depicted in Fig. 6. Then, incorporating the correct behaviour of the potential at 
small distances, we finally obtain the potential energy in the following form 



Certainly, with increasing time the mean distance between particles increases, and the 
Coulomb potential is thus switched on at distances larger than the freeze-out particle mean 

— 1/3 

distance n f . But we restrict our calculation to the freeze-out mean distance. It means 
that in the competition between the two potentials fTcouiomb and U cxpan we overestimate the 
contribution of the distorted potential {7 expan and hence we overestimate the influence of the 
multiparticle environment. 

The results of the calculation with potential ( |37j ) are shown for the LHC freeze-out 
conditions as solid curves in Fig. 4. The same ratio G /G COT for the SPS freeze-out conditions 
is given in Fig. 7. The correction factor we obtained reveals only small deviations from the 
standard Gamov factor. This result can be explained by a fast decrease of the density of 
secondary particles with increasing distance of the pair from the fireball, which in turn results 
in a very small distortion of the two-particle Coulomb potential. 



1 



-tanh [2{x - 2)] , 



(36) 



2 




(37) 
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It is instructive to point out that evaluations that were made by the authors in on the 
basis of a quasi-classical approximation are in good qualitative agreement with the present 
results. 

4 Summary and conclusions 

In our first approach (static scenario) it was assumed that the whole position space is filled 
by secondary pions with constant density. This uniform environment of secondary pions 
results in a screened two-pion Coulomb potential. We showed that for future LHC and 
RHIC experiments the screening radius of the Coulomb interaction at the freeze-out density 
is of a size comparable with the source, and therefore the factorization of Eq. (|]) || is no 
longer valid. Moreover, solutions of the Schrodinger equation produce a correction factor 
G COT = iV'qC 1 ' — 0)| 2 which noticeably deviates from the standard Gamov factor (see Fig. 1). 
However, as we showed further, this model is a quite unrealistic approximation and is not 
relevant to the real picture of an expanding pion system after freeze-out. 

The conclusions reached with the first model change drastically after passaging to a more 
realistic model of an expanding pion system. In the second model, we first reduce the time 
evolution of the multipion post-freeze-out environment to a stationary one, parametrizing 
the density of the secondary pions n n (t,R) for all times t as n(R) = const/i? 2 , where R is 
the distance from the fireball. The 'constant' is determinded by the particular freeze-out 
conditions, namely it should be normalized on the real pion density ri{ at the time of freeze- 
out. This parametrization results from the inequality n n (t, R) < const/i? 2 , where equality is 
reached for an expanding system in which the radial particle velocities are all equal. Quite 
obviously, owing to the particle velocity dispersion, this parametrization is an overestimation 
of the real post-freeze-out pion density at any time t. Consequently, adopting this model 
of the pair environment, we overestimate the deviation from the Coulomb potential and 
thus overestimate the deviation from the Gamov factor. A further reduction results from 
the time evolution of the mean radii: R(t), which is the distance of the separate pion pair 
c.m.s. from the fireball centre, and r(t), which is a classical distance between these pions. 
Exploiting that reduction, we parametrize the relative motion, as well as the pair c.m.s. 
motion, by one variable r, thus eliminating the time dependence. After these reductions 
only one independent variable r is left. 

It should be pointed out that the time derivatives can also be included into our consid- 
eration. From the r-parametrization, we obtain a proportionality of the time derivative to 
the pion-pion relative velocity w re i, namely = v m \ J^. Since we are interested only in the 
small relative pion momenta this derivative can be neglected for a first-order estimate. 

So, we reduce the problem to a stationary one where, in contrast to the first scenario, the 
pion pair moves after freeze-out in a non-uniform environment. Effectively, this means that 
we consider the problem in the two-pion rest frame, where the relative radial motion of the 
particles that create the environment is slow enough with respect to the radial expansion. 
Practically it allows us to consider the stationary Schrodinger equation instead of the time- 
dependent one. We show that the main contribution to the behaviour of the correction factor 
comes from the behaviour of the potential at the large distances that separate interacting 
particles. On the other hand, the interacting particles are separated by large distances when 
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they are far enough from the fireball (this statement is a basis of the r-parametrization). 
Hence, because of a very fast decrease of the density with respect to the distance from the 
fireball R, namely n(R) = const /R 2 , the density of secondaries is small or even negligible 
at these distances. At this stage of the evolution, the long-range part of the potential is 
just the Coulomb one and it is responsible for the behaviour of the correction factor at 
small relative momentum q < 50 MeV/c. That is why the correction factor for this region 
of relative pion momentum practically coincides with the Gamov one. The short-range 
behaviour of the potential (when r is smaller than a mean distance between the particles 
that form the environment), which is the Coulomb one for point-like pions, provides also 
only small deviations of the correction factor G cor from the standard Gamov factor G for a 
large pion relative momentum q >50-70 MeV/c, as can be seen in Figs. 4 and 7. 

For the high "Li^C" freeze-out pion density rif = 0.25fm _1,/3 the reduction of the correc- 
tion factor G cor , as seen in Fig. 4, to the standard Gamov factor G increases with decreasing 
Viei/v cm , the ratio of the relative velocity of the detected pions to their cm. velocity in the 
fireball rest frame. When this parameter is much less than unity, the pion pair promptly 
escapes the initial high-density region so that slow relative motion of the two pions takes 
place in an approximately empty spatial region and the distortion of the mutual Coulomb 
potential is weak. 

This is not the case for "^PiS" freeze-out conditions, as seen in Fig. 7. For the compare- 
tively low pion density rif = 0.03 fm~ 1//3 the effect of small relative velocity is not pronounced 
because from the very beginning of the freeze-out the pions are not in such a dense envi- 
ronment and even small cm. velocities are quite sufficient to bring the pair promptly into 
regions where the influence of the environment is negligible. 
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Figure 1: Correction factor G cor (q) =\ "0 q (r = 0) | 2 as a function of the relative pion 
momentum q = |q| (MeV/c). Wave function ^q(r) is the solution of the Schrddinger equation 
for the screened Coulomb potential (10). Three curves correspond to different freeze-out 
conditions: 1) LHC freeze-out conditions: rif = 0.25 fm -1 / 3 , Tf = 190 MeV, and i? scr = 7.9 
fm (see Eq. (9)); 2) Intermediate case: R scr = 19.3 fm; 3) SPS freeze-out conditions: rif = 
0.03fm~ 1/3 , T f = 190 MeV, and R scv = 22 A fm. 
Bottom curve is the standard Gamov factor, i? scr = oo. 
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Figure 2: Dependence of the spherical pion density n sp h(£,-R) = 4nR 2 n(t, R) (fm -1 ) on R, 
where n(t, R) is the number of pions in the unit volume at time t and at distance R from 
the center of the fireball (J °° dRn sph (t, R) = N w ). The density was evaluated for initial data 
associated with SPS freeze-out conditions: rif = 0.03fnT 1/3 , Rt = 7.1 fm, T f = 190 MeV. 
The spherical pion density (solid Gaussian-like curves) is depicted for the times: 1) t = 0, 
initial distribution, 2) t = 10 fm/c, 3) t = 30 fm/c, 4) t = 70 fm/c, 5) t = 110 fm/c, 6) 
t = 140 fm/c, 7) t = 170 fm/c, 8) t = 200 fm/c. The dotted curve is the evolution of the 
maximum of spatial pion distribution. The horizontal solid line at the top is a constant 
spherical density n sp h(i, R) = AitR 2 n(R) = const (fm -1 ). 
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Figure 3: Exponent b(q) (see Eqs. (31) and (32)) 

a) for the LHC freeze-out conditions: m = 0.25fnT 1/3 , T f = 190 MeV; 

b) for SPS freeze-out conditions: n { = 0.03fnT 1/3 , T f = 190 MeV. 
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Figure 4: The ratio of the standard Gamov factor Gq to the correction factor G COT (q) 



o) 



obtained from the solution of the Schrodinger equation for the LHC freeze-out 
conditions: n { = 0.25fm" 1/3 , R { = 7.1 fm, T f = 190 MeV. 
The dotted curves show the evaluation made using the potential (36), the solid curves using 
the improved potential (38). The curves are drawn for different mean momenta of the pion 
pair P cm = | p a + p&|/2 in the fireball frame. 



15 




10 20 30 40 50 60 70 

Relative pion momentum q (MeV/c) 



Figure 5: Dependence of the dynamical freeze-out radius r(q) = Rf^ 1 (fm) on the relative 
pion momentum q (MeV/c). Slope lines are drawn for different mean momenta of the pair: 

— 1/3 

-Pcm — 50, 150, 450 MeV/c. The mean distance between particles at freezeout a — n { 
(fm) are drawn as horizontal lines for SPS (a = 3.16 fm) and LHC (a = 1.58 fm) freeze-out 
conditions. 
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Figure 6: The solid curve is a smooth potential switcher s(x) — \ — |tanh [2(x — 2)]. The 
dotted curve is an alternative potential switcher 1 — s(x). 
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Figure 7: The ratio of the standard Gamov factor G to the correction factor G cor (q) = 



ipq(r = 0) | 2 obtained from the solution of the Schrodinger equation, using the improved 
for SPS freeze-out conditions: m = 0.03far 1/3 , P f = 7.1 fm, T f = 190 MeV. 
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